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Recent experiments in optimally hole-doped iron arsenides have revealed a novel magnetically 
ordered ground state that preserves tetragonal symmetry, consistent with either a charge-spin density 
wave (CSDW), which displays a non-uniform magnetization, or a spin-vortex crystal (SVC), which 
displays a non-collinear magnetization. Here we show that, similarly to the partial melting of the 
usual stripe antiferromagnet into a nematic phase, either of these phases can also melt in two 
stages. As a result, intermediate paramagnetic phases with vestigial order appears: a checkerboard 
charge density-wave for the CSDW ground state, characterized by an Ising-like order parameter, 
and a remarkable spin-vorticity density-wave for the SVC ground state - a triplet d-density wave 
characterized by a vector chiral order parameter. We propose experimentally detectable signatures of 
these phases, show that their fluctuations can enhance the superconducting transition temperature, 
and discuss their relevance to other correlated materials. 


I. INTRODUCTION 

One of the hallmarks of the superconducting state of 
the iron-based materials [l| is its typical proximity to 
a stripe magnetically ordered state, with spins aligned 
parallel to each other along one in-plane direction and 
anti-parallel along the other (see Fig. [T^) 0. As a re¬ 
sult, this stripe state breaks two distinct symmetries of 
the high-temperature paramagnetic-tetragonal state: a 
continuous spin-rotational 0(3) symmetry and an Ising- 
like Z 2 symmet^ related to the equivalence of the x and 
y directions [J-Q- Magnetic fluctuations present in the 
paramagnetic state can cause these two symmetries to be 
broken at different temperatures, giving rise to an inter¬ 
mediate nematic phase that preserves the spin-rotational 
0(3) symmetry but, as a “vestige” of the stripe order 
[T^ . breaks the tetragonal Z 2 symmetry [Tl|. Indeed, in 
the phase diagrams of most iron-based superconductors, 
the magnetic transition line is closely followed by the 
structural/nematic one at slightly higher temperatures. 
The corresponding nematic degrees of freedom impact 
not only the normal state electronic properties (12 h^ 
but also the onset and gap structure of the supercon¬ 
ducting state (22l - [^ . 

Recently, experiments in the hole-doped pnictides 
Ba(Fei_a;Mna;) 2 As 2 (Bai_a:Na3;)Fe2As2 ( 2 ^, and 

(Bai_xKa;)Fe 2 As 2 have revealed another type of 
magnetically ordered state that does not break the 
tetragonal Z 2 symmetry of the lattice. Neutron scatter¬ 
ing experiments showed that its magnetic Bragg 

peaks are at the same momenta as in the stripe magnetic 
phase - namely, Qi = (tt, 0) and Q 2 = (0,7r) in the Fe- 
only Brillouin zone. Consequently, it has been proposed 
[28 - [^ that the tetragonal magnetic state is the re¬ 
alization of one of two possible biaxial (i.e. double-Q) 
magnetic orders [^, l3^ - [^ . One possibility is a “charge- 
spin density wave” (CSDW), displaying a non-uniform 
magnetization which vanishes at the even lattice sites 


and is staggered along the odd lattice sites (Fig. [Ud). 
The other option is a “spin-vortex crystal” (SVC), in 
which the magnetization is non-collinear (but coplanar) 
and forms spin vortices staggered across the plaquettes 
(Fig. [T}:). Both CSDW and SVC phases are tetragonal, 
but have a unit cell four times larger than the param¬ 
agnetic phase. Interestingly, in (Bai_a,Naa,)Fe 2 As 2 and 
(Bai_a,Ka;)Fe 2 As 2 , the tetragonal mag netic state is ob¬ 
served very close to optimal doping (^ . 1^ , where super¬ 
conductivity displays its highest transition temperature. 
Therefore, understanding the properties of these biaxial 
tetragonal magnetic phases is important to assess their 
relevance for the superconductivity. 

In this paper, we show that both the CSDW and the 
SVC magnetic phases support composite order parame¬ 
ters that can condense at temperatures above the onset 
of magnetic order, and whose fluctuations can help en¬ 
hancing Tc- As with the nematic phase, these partially 
ordered phases are paramagnetic, i.e. fluctuations re¬ 
store the time-reversal symmetry that is broken in the 
ground state. In contrast to the nematic phase, how¬ 
ever, they preserve the point group symmetry of the lat¬ 
tice, but break other symmetries, including translational 
symmetry (^ . In particular, upon melting the CSDW 
phase, we find a vestigial Ising-like charge-density wave 
(CDW) phase with ordering vector Qi -I- Q2 = (tJ’iTt), 
in which the previously magnetized sites acquire a dif¬ 
ferent charge than the previously non-magnetized sites. 
On the other hand, upon melting the SVC ground state, 
we find a vestigial phase that retains memory of the pre¬ 
ferred plane of magnetization (in spin space), and of the 
staggering of the spin vortices across the plaquettes. This 
spin-vorticity density-wave (SVDW) is a triplet d-density 
wave characterized by a vector chiral order parameter, 
which is manifested as a spin-current density-wave with 
modulation Qi + Q2 = (tTjTt). Besides shedding light on 
the magnetism of hole-doped iron pnictides, our results 
provide a novel microscopic mechanism for the forma¬ 
tion of d-density waves, which have also been proposed 
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Figure 1: Magnetic ground states of the iron pnictides: 
(a) stripe antiferromagnet, (b) charge-spin density-wave 
(CSDW), (c) spin-vortex crystal (SVC). The first is or¬ 
thorhombic with a doubled unit cell; the latter two remain 
tetragonal but with a quadrupled unit cell. Mi and M2 are 
the magnetic order parameters corresponding to the ordering 
vectors Qi = (7r,0) and Q2 = (0, tt). The left panels are 
the actual spin density-wave patterns in real space, whereas 
the right panels are schematic representations focusing on the 
magnetization at the lattice sites. 


in cuprates [s^ and heavy fermions M- 

The paper is organized as follows: in Section |ll] we 
present the theoretical model that gives rise to the ves¬ 
tigial CDW and SVDW orders. Section Hill discusses the 
implications of these vestigial orders for both the normal 
state and superconducting state properties. Concluding 
remarks are presented in Section To make the pa¬ 
per transparent and accessible, all formal details are pre¬ 
sented in appendices. Appendix contains the deriva¬ 
tion of the saddle-point equations that give the phase 
diagram of the SVDW phase discussed in Section HIl In 
Appendix |B] we derive microscopically the free energy dis¬ 
cussed in Section IIIII Finally, Appendix [0 presents the 
derivation of the effective pairing interactions promoted 
by CDW and SVDW fluctuations discussed in Section 

Em 


II. THEORETICAL MODEL FOR THE 
VESTIGIAL PHASES 

A. Effective action 

We define two magnetic order parameters. Mi and 
M 2 , associated with the two ordering vectors Qi = (tt, 0) 
and Q 2 = (0,7r), respectively. Thus, the local spin is 
given by S (r) = Mie®*^' '’. As discussed in Refs, 

[a, H, [28l - [^ [32 - [^ , the most general lowest order action 
that respects the tetragonal and spin-rotational symme¬ 
tries is given by: 

5 [M,] = ^ X,-' {Ml + M 2 ) + I ^ {Ml + Ml )" 

" I / (Ml ■ ^ 2 )" + (Ml ■ M 2 )'' ■ (1) 

For simplicity, we will consider the finite temperature 
problem, but the same conclusions can be extended to 
the quantum case. Here, f^ = f and f^ = f d^x 
where q is the momentum and x is the position. In the 
neighborhood of a finite T magnetic transition, and for 
a quasi-2D system, we can use the small q expansion 
Xq^ ~ ro+ql\+Jz sin^ where xq is the distance to the 
mean-field magnetic critical point. 

The quartic coefficients u, g, w determine the nature 
of the magnetic ground state. These are, in turn, sensi¬ 
tive to microscopic considerations. The localized J 1 -J 2 
model favors positive g and w [131. On the other hand, 
itinerant approaches (at weak and strong coupling) have 
found parameter regimes in which q and w can be ei¬ 
ther positive or negative H, [Isl - lloL H^. For 

g > max (0, — w), the energy is minimized by the stripe 
state shown in Fig. [1^, in which either (Mi) = 0 or 
(M 2 ) = 0. Thus, in addition to breaking the 0(3) spin- 
rotational symmetry, the magnetic ground state sponta¬ 
neously breaks a Z 2 symmetry by selecting one of the two 
order parameters to be non-zero. Since Mi and M 2 are 
related by a 90° rotation, once this Z 2 symmetry is bro¬ 
ken the tetragonal symmetry of the system is lowered to 
orthorhombic (see Fig. A composite Ising-nematic 
order parameter, living on the bonds of the lattice, can be 
identified by performing a Hubbard-Stratonovich trans¬ 
formation on the quartic term with coefficient g, yielding 
(</?nem) = <7 (Mi — M|). Because Z 2 is a discrete symme¬ 
try, while spin-rotational 0(3) is a continuous symmetry, 
a strongly anisotropic 3D system will generically display 
a vestigial paramagnetic nematic phase where (MJ = 0 
but ((^nem) iim. 

For g < max(0, —w), the ground state of Eq. ([T|) is 
no longer a uniaxial magnetic stripe state, but a biax¬ 
ial magnetic state with | (Mi) | = | (M 2 ) | that preserve 
tetragonal symmetry. If tc < 0, the energy is minimized 
by (Ml) II (M 2 ), which in terms of the local spin config¬ 
uration S (r) corresponds to a non-uniform state as de¬ 
picted in Fig. [Dd. We identify this state as a charge-spin 
density-wave (CSDW). On the other hand, if w > 0, the 
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Figure 2: The vestigial composite states associated with (a) 
the stripe antiferromaget, (b) the CSDW state, and (c) the 
SVC state. The real-space spins (in gray) and the magnetic 
order parameters in spin space (red and blue arrows, repre¬ 
senting Ml and M 2 , respectively) should be understood as 
fluctuating, i.e. (Mi) = 0 in all cases. In (a), the vestigial 
state is nematic (unequal blue and red bonds), associated with 
selecting between Mi and M2 fluctuations in spin space. The 
original unit cell is shown as a dashed square. In (b), the ves¬ 
tigial state breaks translational symmetry via a checkerboard 
charge density-wave (unequal blue and red sites). Mi and M2 
are locked to be collinear in spin space. In (c), the vestigial 
spin-vorticity density-wave state breaks inversion and trans¬ 
lational symmetries via a staggered pattern of spin vortices 
in the center of the plaquettes (unequal blue and red plaque- 
ttes). The corresponding spin-current pattern is shown by the 
green arrows. Because Mi and M2 are locked to be orthog¬ 
onal in spin space, the residual spin-rotational symmetry is 
0(2) instead of 0(3). Both (b) and (c) preserve tetragonal 
symmetry, as shown by the dashed-line unit cell. 


energy minimization gives (Mi) _L (M2), correspond¬ 
ing to a non-collinear, coplanar spin configuration (see 
Fig. [It). This state is identified as a spin vortex-crystal 
(SVC). We now discuss whether these tetragonal mag¬ 
netic phases can melt in a two-stage process, giving rise 
to vestigial orders akin to the nematic phase. 


B. Charge-spin density-wave 

Consider the CSDW state: Once the magnetization 
direction is chosen by spontaneous breaking of the 0(3) 
spin-rotational symmetry, there remains a four-fold de¬ 
generacy corresponding to whether Mi and M 2 are par¬ 
allel or anti-parallel to the chosen direction. As is ap¬ 
parent in Fig. [Dd, this corresponds to the breaking of 
translational symmetry, leading to a four-site unit cell. 
Notice, however, that the product of a translation by 
the vector x -|- y followed by time-reversal is preserved. 
Thus, there is an essential Z 2 symmetry that interchanges 
the magnetic and non-magnetic sublattices of the CSDW 
state. 

The order parameter field for this Z 2 symmetry is ob¬ 
tained via a Hubbard-Stratonovich transformation on the 
quartic term with coefficient 2w in Eq. U]), ((/3 cdw) = 
2 |w| (Ml • M 2 ). Clearly, :/3 cdw is a scalar that carries 
momentum Qi + Q 2 = (tt, tt), i.e. the condensed phase 
is a CDW that doubles the unit cell, but leaves time- 
reversal and the tetragonal symmetry of the lattice intact 
(see Fig. j^h)- Thus, in real space, the CDW order pa¬ 
rameter lives on the lattice sites. The fact that the unit 
cell decreases from four to two sites upon going from the 
CSDW to the CDW phase is due to the restoration of 
time-reversal symmetry, which implies the restoration of 
the translational symmetry by x -|- y. A simple change 
of variables in Eq. HI), Ml ^ 2 - 1/2 [Ml -b M2] and 
Ml ^ 2 -i/2[Mi-M 2], interchanges the identities of the 
two scalar orders, :/3nem ^ <PCDW, but leaves the form of 
S unchanged albeit with {g,w) —>■ — (w,g). Thus, the 
properties of the CDW phase are akin to those of the 
Ising-nematic phase - in particular, a quasi-2D system 
will again display for a range of intermediate tempera¬ 
tures a phase with (M^) = 0 but ((/^CDw) 7^ 0. 


C. Spin-vorticity density-wave 

Consider now the SVC state, characterized by two 
equal magnitude orthogonal vectors Mi and M 2 . Upon 
fixing the direction of Mi, which breaks the 0(3) spin- 
rotational symmetry, there remains an additional 0(2) 
symmetry related to choosing M 2 in any direction along 
the plane perpendicular to Mi [i^. Thus, the SVC phase 
can be completely characterized by a pseudo-vector or¬ 
der parameter V^svdw that specifies the ordering plane 
which contains Mi and M 2 , and also by the orienta¬ 
tion of Ml within that plane. <PsvDW obtained via 
a Hubbard-Stratonovich transformation of the quartic 
term w (Mi • M 2 )^ —>■ —w (Mi x M 2 )^ in Eq. H]), yield¬ 
ing (<PsvDw) = (Ml X M 2 ), which can be identified 
as a vector chirality [HI, |4^. Thus, upon approaching 
the SVC phase from high-temperatures or by melting it, 
there can be an intermediate state where (^’svdw) 0 
but the orientation of Mi is not fixed, (Mi) = 0. This 
chiral paramagnetic state preserves time-reversal sym¬ 
metry and retains the memory of the staggering pat- 
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Figure 3: Phase diagram, within the saddle-point approxi¬ 
mation, of the coupled SVDW paramagnetic and SVC mag¬ 
netic transitions. In the shaded area, where the out-of-plane 
anisotropy is strong, the two transitions are split. The tun¬ 
ing parameters are the Ginzburg-Landau coefficients u and 
w (see Eq. [TJ and the magnetic coupling between layers Jz. 
The parameter w is given by ui = Tn,qw/2tt, as discussed in 
Appendix IbI 


tern of spin vortices along the plaquettes in the SVC 
phase, and is therefore called a spin-vorticity density- 
wave (SVDW) [12. Note that the vector chiral order 
parameter produces an emergent Dzyaloshinskii-Moriya 
coupling (PsvDW ■ (^1 ^ ^ 2 ) relating the translational 
symmetry-breaking to a preferred “handedness” in spin- 
space. In the SVDW state, not only is the translational 
symmetry lowered by the doubling of the unit cell (since 
V’svDW carries momentum Qi -|- Q 2 = (7r,7r)), but also 
the soft spin fluctuations near the magnetic transition are 
constrained to lie in the plane defined by <PsvDW 

m- 

Because Vsvdw breaks a continuous 0(3) symmetry, 
there are two Goldstone modes in the SVDW phase. 
Consequently, in contrast to the Ising-nematic cases, the 
Mermin-Wagner theorem does not ensure the existence 
of the SVDW phase even in the two-dimensional limit. 
To investigate whether (<PsvDw) ^ b while (Mi) = 0 
is possible, we calculated the phase diagram for a mag¬ 
netic SVC ground state treating the action in Eq. © in 
the saddle point approximation (see Appendix E|. We 
find that for a strongly anisotropic system, i.e. Jz -C w, 
there is a wide range of values oiu/w for which there are 
two transitions, with an intermediate SVDW phase and 
a low-temperature SVC phase (see Fig. [3]). However, in 
this approximation, the transition to the SVDW phase is 
always first-order. 

Spin rotational symmetry is not an exact symmetry of 
nature, and indeed most iron pnictides display a sizable 
spin anisotropy (43l - [45l |. Because the ordered moments 
tend to point parallel to the FeAs plane, the most signif¬ 
icant effects of spin-orbit coupling can be captured phe¬ 
nomenologically in Eq. m by including an easy-plane 
anisotropy term k (M^ ^ -|- z) with coupling constant 
K > 0 |46l | . The spin rotational symmetry is thus reduced 


to 0(2) and the SVDW chiral order parameter becomes 
the pseudo-scalar (/Jsvdw = 2w (Mi x M 2 )-z, which only 
breaks a discrete chiral Z 2 symmetry. For such 0(2) x Z 2 
model, it is known from both numerical and analytical 
investigations that in 2D the Z 2 symmetry is broken at 
higher temperatures than the Kosterlitz-Thouless transi¬ 
tion of the 0 ( 2 ) order parameter [13,13, there is 
no doubt that there is a vestigial chiral SVDW phase. 
The extent to which the spin anisotropy is quantitatively 
significant depends on the (currently unknown) value of 
the ratio k/ (Tsvdw — Tsvc)- 


III. MICROSCOPIC IMPLICATIONS OF THE 
VESTIGIAL ORDERS 


A. Normal-state manifestations 

To discuss the experimental manifestations of the ves¬ 
tigial CDW and SVDW states, we investigate their cou¬ 
pling to the low-energy electronic states of the pnic¬ 
tides. We consider a three-band model (^, § with 
a circular hole pocket k at the center of the Bril- 
louin zone, and two elliptical electron pockets ■fei.k-eQi 
and Ce 2 .k+Q 2 centered at momenta Qi = (tt, 0) and 
Q 2 = ( 0 , tt), respectively (see Fig. il). The magnetic 
order parameters couple to these electronic states via 

Z(k,a /3 ■ (^ap (4.kQCe,,k/3 + h-C.) , where the operator 

Ca.ka annihilates an electron in band a with momentum 
k (measured with respect to the center of the pocket) 
and spin a, and (Tap are Pauli matrices. We further in¬ 
troduce magnetic A 5 and charge Ac order parameters 
with ordering vector Qi -|- Q 2 = ( 7 r, 7 r), which couple to 
the electronic states via 


■ 0-«/jc);^,kaCei,k/3 + ^-C.] , 

k,Q:/9 

nc=J2 [^Ci5„^cj;^,kaCe,,k/3 + h-C. 

k,a^ 


( 2 ) 


Here these fields have real and imaginary parts, Ag = 
Ag + zAg and Ac = AJ^ -I- iA^, where the real parts 
correspond to conventional SDW or CDW orders, while 
the imaginary parts corresponds to spin or charge cur¬ 
rent orders. By integrating out the electronic degrees of 
freedom, we obtain the coupling between M^ and Ag, 
Ac to lowest-order in the action (see Appendix |B|) : 


,55eff = A [A'^ • (Ml X M 2 ) - A^ (Ml • M 2 )] (3) 

with the coefficient A = i J^Gh,kGej,kGe2,kj where 
^a \ ~ ~ Corresponding non-interacting 

Green’s function. As expected, the Ising-like order pa¬ 
rameter (/5CDW oc Ml • M 2 induces a checkerboard-like 
charge order (see Fig. [ 3 } 3 ). On the other hand, the 
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Table I: Magnetic ground states of the pnictides and their corresponding vestigial states. Mi and M 2 are the magnetic order 
parameters corresponding to the ordering vectors Qi = (tt, 0) and Q2 = (0, tt). 


magnetic ground state 

vestigial state 

broken symmetry 

real space pattern 

physical manifestation 

stripe: (M 2 ) or (Mi) = 0 

nematic: (Mf — M^) 7 ^ 0 

rotational (tetragonal) 

unequal bonds 

orthorhombic distortion 

CSDW: (Ml) II (M 2 ) 

CDW: (Ml ■ M 2 ) / 0 

translational 

unequal sites 

charge density-wave 

SVC: (Ml) T (M 2 ) 

SVDW: (Ml X M 2 ) / 0 

translational + inversion 

unequal plaquettes 

spin-current density-wave 



Figure 4: Schematic Fermi surface of the iron pnictides, with a 
central hole pocket and elliptical electron pockets. The wavy 
lines represent the inter-pocket pairing interactions generated 
by the magnetic fluctuations (repulsive > 0 ) and by fluc¬ 
tuations of the vestigial CDW and SVDW states (attractive 
U < 0). 


SVDW order parameter 93 svdw ^1 ^ ^2 is mani¬ 
fested as a spin-current density-wave with propagation 
vector (tTjTt), i.e. a spin current polarized parallel to 
V’svDW propagating along the bonds of the lattice 
in a staggered pattern across the square plaquettes (see 
Fig. [ 2 j;). Thus, the SVDW corresponds to a triplet d- 
density wave |^. 

Note that probing the CDW via x-rays may be difficult, 
since the hybridization between Fe and As/Se doubles 
the unit cell of the Fe-only square lattice, making (tt, tt) 
a lattice Bragg peak. While the real CDW could in prin¬ 
ciple be detected experimentally by a probe sensitive to 
the local charge on the Fe sites, such as STM, detecting 
a spin-current density-wave would be rather challenging. 
Alternatively, one can consider the effects of a Zeeman 
field H. Despite not coupling to (pnem, we find that it cou¬ 
ples to both (/ 5 CDW S'lid V’svDW the action via the terms 
yA^H • (Ml X M2) and 7 (H • Ag) (Mi • M2), with the 
same Ginzburg-Landau coefficient 7. Therefore, in the 
presence of a magnetic field, a pattern of staggering or¬ 
bital currents (i.e. a singlet d-density wave | 36 l |l appears 
in the SVDW state, which can in principle be detected 
by NMR. Table U summarizes the magnetic ground states 
of the pnictides along their vestigial paramagnetic states. 


B. Impact on the superconducting state 

Fluctuations of the SVDW and CDW states arise from 
four-spin correlations, and are complementary to the 
magnetic fluctuations that arise from two-spin correla¬ 
tions. An important issue is whether these fluctuation 
modes promote compatible superconducting states. Be¬ 
cause the magnetic fluctuations are peaked at momenta 
Qi = (tt, 0 ) and Q2 = (0,7r), they promote a repulsive 
inter-pocket interaction V > 0 between the hole and 
the electron pockets (see Fig. S]). Solution of the cor¬ 
responding linearized gap equations yields the so-called 
s"' state, where the gap functions have different signs 
in the electron and in the hole pockets [i^. The tran¬ 
sition temperature is given by Tc oc exp^—with 

Aq = y/ 2 NhNf, V, and Na denoting the density of states 
of band a. 

The SVDW and CDW fluctuations, on the other hand, 
are peaked at the momentum Qi -|- Q2 = (7r,7r) and 
promote an attractive inter-pocket interaction U < 0 
between the two electron pockets (see Fig. | 4 ] and Ap¬ 
pendix P. Solution of the linearized gap equation re¬ 
veals that the leading eigenstate remai ns the s"' one, 
but the eigenvalue is enhanced, A = + Ajj, 

with Ajj = AAV > Q Therefore, fluctuations associated 
with these vestigial states may enhance the value of 
promoted by spin-fluctuations pairing, without affecting 
the symmetry of the Cooper pair wave-function. Similar 
conclusions have been found for the combination of pair¬ 
ing promoted by nematic fluctuations (peaked at Q = 0 ) 
and magnetic fluctuations (peaked at Qi and Q2) [^ . 

IV. CONCLUDING REMARKS 

In summary, we showed that both biaxial tetrago¬ 
nal magnetic ground states of the pnictides - the non- 
uniform CSDW and non-collinear SVC states - can melt 
in two-stage processes, giving rise to CDW and SVDW 
vestigial states, respectively. While both preserve the 
point-group and time-reversal symmetries, but break the 
translational symmetry of the iron square lattice, only 
the SVDW state also breaks inversion symmetry by en¬ 
tangling the spin-space handedness to a doubling of the 
real-space unit cell. Because in the iron superconductors 
the hybridization with the puckered As atoms already 
doubles the unit cell of the Fe square lattice, the CDW 
and SVDW states are more rigorously classified as intra- 
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unit-cell orders. Recent experiments on Sri_xNaxFe 2 As 2 
(S^ l and Bai_xKxFe 2 As 2 [l^ found direct evidence for 
a low-temperature CSDW phase, which can support a 
CDW vestigial phase. It remains to be seen whether the 
tetragonal magnetic phase can be reached in these com¬ 
pounds without first crossing the stripe magnetic state. 
In contrast, in Ba(Fei_a;Mn 3 ,) 2 As 2 [^, the tetragonal 
magnetic phase has been reported to exist over a wide 
doping range as the primary instability of the paramag¬ 
netic phase. 

Beyond the physics of iron-based superconductors, our 
results establish the melting of double-Q orders as a mi¬ 
croscopic mechanism to create d-density wave states. The 
latter have been proposed to be realized in other strongly 
correlated systems, such as the pseudogap phase of un¬ 
derdoped cuprates and the hidden-order phase of 
the heavy fermion compound URu 2 Si 2 [sol - f^ . mostly 
on phenomenological basis. Whether our mechanism is 


directly applicable to those systems is an appealing topic 
for future investigation. 
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yielding: 


5 [M,] = ^ {Ml + M|) + ^ {Ml + M|)' - ^ {Ml - M|)' - 2n; £ (Ml X M 2 )' (A3) 

Hereafter for simplicity we introduce the parameters g = g + w and u = u + w. Since we are interested in the 
vestigial phase of the spin vortex-crystal, which has tetragonal symmetry, the nematic order parameter (M^) — (M|) 
never condenses, and we can ignore the corresponding quartic term. Introducing the Hubbard-Stratonovich fields 
corresponding to the other two quadratic terms, we obtain: 

g2u.(MixM2) (MixM2) = M J (A4) 

Here, V’svdw spin-vorticity density-wave (SVDW) vectorial order parameter whose mean value is given by 

(tpgVDw) = 2i/; (Ml X M 2 ). The field tp is not an order parameter, and just renormalizes the magnetic correlation 






length via {ijj) = u i.e. it corresponds to Gaussian magnetic fluctuations. Thus, the effective action 

is given by: 


5 [M„ <PsvDw] = ^ (x,-' + V') {m! + Mi) - 2 £ C^SVDW ■ (Mi x M 2 ) + ^ (A5) 

Approaching the SVDW phase from the paramagnetic state, we can integrate out the magnetic degrees of freedom, 
yielding an effective action for ip and y’svDW- 


SeS [4’t VsvDw] — 


‘/^ivDW , 1 


210 




IA6) 


where are the eigenvalues of the matrix Ay corresponding to the Gaussian action in M^. The Gaussian part of 
the action can be rewritten in the convenient matrix form: 


(Mi M 2 ) 


Xq^ +4 

0 

0 

0 

-Pz 

Py ^ 

0 

x4^ +4 

0 

Pz 

0 

Px 

0 

0 

X^^ +'0 

-Pv 

Px 

0 

0 

•■Pz 

-Pv 

x4^ +4 

0 

0 

-<^z 

0 

Px 

0 

xli^ + 4 

0 


-•Px 

0 

0 

0 

x4^ + 4 ) 

gives: 







Ml 

M 2 


(A7) 


SeS [tp, ¥’SVDW] — 


‘T’ivDW 


2w 


2u 


+ / log [{Xq ^ + 'fp) {Xq^ + fp + ‘PSVDw) (Xg ^ ^ - <PSVDw)] 

J a 


(A8) 


So far our result is exact. To proceed, we employ the saddle-point approximation to determine the equations of 
state for %p and ^/JsvdWi which corresponds to self-consistently accounting for the Gaussian magnetic fluctuations. 
The saddle-point equations become: 


‘/’SVDW 


1 


w 


ia Xq 


1 


t. 

u 


Xq^ +4’ - ‘PSVDW Jq Xq^ + 4 + ‘/’SVDW 

1 f 1 

—-^ - ' / —1 

Xq +4 - ‘/’SVDW jq Xq +4 + ‘/’SVDW Jq Xq +4 


+ 


I 


(A9) 


Since our focus is on the proximity to a finite-temperature magnetic transition, we ignore the spin dynamics and 
use the low-energy expansion for the spin susceptibility appropriate for anisotropic layered systems: 


Xq^ =ro+q'\\+Jz sin^ ^ 


(AlO) 


where ro = a (T — Tjv), a > 0, T/v is the mean-field magnetic transition temperature, = qi + qy, and Jz is the 
inter-layer magnetic coupling. Defining the renormalized magnetic mass: 


r = ro -f V” ‘X C 

where ^ is the magnetic correlation length, we obtain: 


(All) 


‘/’SVDW = W 


\f 1 

f 1 j 

Jq r + q^^ + Jz sin^ ^ - /^svDW 

Jq r + qf^ + Jz sin^ ^ -f (/JsvDW 


(A12) 


r = rQ + u 


1 1 , 

r 1 1 

r 1 1 

Jq ''' + ^1 + Jz sin^ ^ - ‘/’SVDW ^ 

qr + q^^ + Jz sin^ ^ -f /^SVDW ^ 

qr + qf^+Jz sin^ ^ 
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The integrals can be evaluated in a straightforward way (we consider only the ujn = 0 contribution to the sum over 
Matsubara frequencies, since we are interested in the finite temperature transition): 


9,1 


-h sin^ ^ a, 


In f ^ f — 

47r Jo 27r J sin 2 _(_„ x 

^ 

4tt Jq 2tt 


' Jz sin^ 

A2 


Jz sin^ ^ + a 


In ' 

27r . 


ln2A — In ^\/Jz + a + y/aJj 
Defining the renormalized critical temperature fg = a(T — T/v) via: 


(A13) 


SmTjv 2A 
ro = To-\ -^— In ■ 


Stt vX 


(A14) 


we obtain the self-consistent equations: 


‘/^SVDW 

r 


wTn , \/Jz+r + ^svDW + \/r + (fisvuw 

-In — - - 

Stt y/Jz+r — :psvDW + y/f — <PSVDW 


(A15) 


uTm 


{y/Jz + r + V^SVDW + y/r + V^SVDw) {y/Jz + r — :psvDW + y/f — <PSVDw) (vXXr -|- y/r) 


For simplicity, we define the renormalized parameters (w, u) = (w, u) ^ as well a.sa=^ = ^ + l and Jz = Jz/w. 
Then the equations can be written as: 


</^SVDW 


= In 


Jz + r + (PSVBW + y/r + (PSVDW 
Jz + r — (fisvBW + y/r — fswBW 

^\/ Jz + r + (PSVDW + y/r + 


(A16) 


= fg — a In 


Jz + r — V^SVDW + y/r — V^SVDW j J-y/r^ 


J1 


"i/2 


where r, fg, and :/3 svdw were rescaled by w as well. The SVDW transition temperature can be obtained by linearizing 
the equations around v^svdw = 0- From the first equation, we obtain the correlation length ri at the SVDW transition: 


FJ2+4-J, 

ri = — - 

2 

which, when substituted in the second equation, gives the SVDW transition temperature rc,svDw: 


(A17) 


?'0,SVDW 


J? + 4 - J, 


3a In 


+ 4 -|- Jz -|- Y -|- 4 — J; 

VITz 


(A18) 


The magnetic transition temperature rg^mag is signaled by the vanishing of the renormalized magnetic mass, i.e. 
the lowest eigenvalue of the Eq. (RTI) . r — <PsvDW- Therefore, it takes place when r reaches the value r 2 determined 
implicitly by: 
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r 2 = In 


\/t/z + 2 r 2 + Y^ 2 r 2 


The magnetic transition temperature is therefore given by: 


(A19) 


ro,mag = r 2 (1 + a) + a In 


V Jz + r 2 + 

'/j'z 


(A20) 


The SVDW and magnetic transitions are split when ro,SVDW > ^o,mag- The region in the J^j parameter space 
where this condition is satisfied corresponds to the shaded area of Fig. 3 in the main text (recall that — = a — 1). 

To determine the character of the SVDW transition, we can expand fo for small :/7 svdw- Substituting r = ri + 
at/JgvDw equation of (IA16I) and expanding for small (/^svdw gives the coefficient of the quadratic term: 


8 + 3 J 2 

“ =-7^ 

12 Y + 4 

Substituting it in the second equation of (IA16I) and collecting the quadratic terms in v^svdw yields: 


(A21) 


ro ivsvBw) ~ 1^0,SVDW 


16 + 3 (2 + a) 

24/^74 


2 

‘/^SVDW 


(A22) 


Therefore, because the coefficient is always positive, the solution with :/5 svdw 0 is achieved at a larger temperature 
than the solution with (/Jsvdw = 0) in other words, fg (</?svdw > 0) > fg (v^svdw — t 0). As a result, the SVDW 
transition is first-order within the saddle-point approximation, even when it is split from the magnetic transition. 


Appendix B: Derivation of the Ginzburg-Landan free energy 


Our starting point is a 3-band model with a circular hole pocket h centered at (0,0) and two elliptical electron 
pockets ei _2 centered at Qi = (7r,0) and Q 2 = (0,7r), respectively. The band dispersions can be conveniently 
parametrized by Q: 


2m 

^ei.k+Qi = ^k — (^0 + ^2 COS20) 

^ea.k+Qa = ■Ck - (^0 - ^2 COS 20) (Bl) 

Here, (5g is proportional to the chemical potential and 62 to the ellipticity of the electron pockets. The angle 9 is 
measured relative to the axis. The non-interacting Hamiltonian is therefore given by (hereafter sums over repeated 
spin indices are implicitly assumed): 


-Ho — ■C/t.kcjj 'Cei,kc7kcr''ei,ko- +^€2 ^ka-^e2 ,k(T (B2) 

k k k 

These electronic states couple to the magnetic order parameters Mi and M 2 according to: 

-Hmag — ^ ^ M^i ■ j^^^TcK/3C/j.k/3 " 1 ” h.C.^ (14^) 

k ,2 

In principle, this last term can be obtained via a Hubbard-Stratonovich transformation of the original interaction 
terms projected into the magnetic channel, as shown in Ref. Q. Here, because we are interested in the higher-order 
couplings of the action involving the M^ order parameters, we neglect these interaction terms, since they only affect 
the quadratic terms of the action. 
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1 . Absence of magnetic field 

In the case where there is no external magnetic field, we focus on the two types of fermionic order that couple 
directly to the SVDW order parameter, Mi x M 2 , and to the CDW order parameter Mi • M 2 . Thus, we introduce 
the Qi + Q 2 = (tt, tt) spin-current density-wave Ag (i.e. a purely imaginary SDW) and the checkerboard charge order 
A'q (i.e. a purely real CDW) defined by: 


'^iS — i ^ ' ®'a/3 (cL,ka'^ei,k/3 .ka'^ea.k/?) 

k 

'^C=J2 (cl2.kaCe,,k/3 + .ka^e, .k/s) (B4) 

k 

To proceed, we introduce the 6 -dimensional Nambu operator: 

4'k = ( cji^kt '^ei.kt '^L.kf '^L.kf *^12.kf ) 

which allows us to write the fermionic action in the compact form: 

S = - [ '4'k + 5o [M^] (B6) 

In the previous expression. So [Mf] corresponds to the terms Mf that arise from the decoupling of the fermionic 
interactions. As we explained above, these terms can be ignored for our purposes. The total Green’s function is given 
by: 


The bare part is: 


e.-' = (sr)‘ 


— Dmag — ViS — Vc 


(B7) 


qT = 


( Gh,k 
0 
0 
0 
0 
0 


V 


0 

Gh,k 

0 

0 

0 

0 


0 

0 

Gei,k 

0 

0 

0 


0 

0 

0 

Gei,k 

0 

0 


0 

0 

0 

0 

Ge2,k 

0 


0 

0 

0 

0 

0 

Ge 2 ,k 


(B 8 ) 


where G^l = iwn — ^i.k are the non-interacting single-particle Green’s functions. The interacting parts are: 


/ 

0 

0 

—Mi^z 

—Mi^x + iMi^y 

— M2,z 

~^2,x + 'i^2,y 

\ 


0 

0 

— — i^l,y 

Mi^z 

— M2,x — *Af2,y 

Af2,2 



—Mi^z 

~Ml,X + 

0 

0 

0 

0 




Mi^z 

0 

0 

0 

0 



— Af2,z 

~M2^X + iM2^y 

0 

0 

0 

0 


V- 

-M2^x — 

M2,z 

0 

0 

0 

0 

/ 


(B9) 

and: 




/ 0 0 0 0 

0 0 0 0 

0 0 0 0 

0 0 0 0 

0 0 —i {Ag^^ — iA'g y) 

y 0 0 -l {A'g^^ + i^s,y) ^^S,z 


0 

0 


* {^S,x + 


0 

0 


0 \ 
0 






0 

0 / 


(BIO) 
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as well as: 


/O 0 0 0 0 0 \ 

0 0 0 0 0 0 

^ _ 0 0 0 0 0 

*^“ 0 0 0 0 0 -A',;, 

0 0 -A'c 0 0 0 

\ 0 0 0 -A^;, 0 0 / 

It is now straightforward to integrate out the fermions and obtain the effective magnetic action: 


(Bll) 


5eff [Mi,M2,A^,Ay 


-Trln 


1 - 






mag 




fc)] »!:>[»»( 





(B12) 


where, in the last step, we expanded for small Mi, M 2 . Here, Tr (• • •) refers to sum over momentum, frequency and 
Nambu indices. A straightforward evaluation gives, to leading order in the coupling between Ag, A^, and M.^: 


5eff [Ml, M 2 , A'^ A'c] = 5 [Ml, M 2 ] + A A" • (Mi x M 2 ) - AA',. (Mi • M 2 ) 


(B13) 


with the coefficient: 


A = 4 


[ Gh,^ 

Jk 


kG 


(B14) 


For perfect nesting. So = S 2 = 0, this coefficient vanishes. For a system in proximity to a finite temperature phase 
transition, expansion in powers of Jq gives: 


A « ApfTY f - - - 2 

+ 6 (fcc„ - e + ^o) 

A « -SSoPfTY f dC—^—— - - - 1 


where pF is the density of states at the Fermi level. Therefore, it is clear that a spin-current density-wave Ag parallel 
to yJsvDW triggered by the SVDW order parameter, <PsvDW ^1 ^ ^ 2 ) whereas a checkerboard charge order A^ 
is triggered by the CDW order parameter :^cdw oc Mi • M 2 . 


7C(3)Pf 

2 f^T 


(B15) 



2 . Non-zero magnetic field 

In the presence of a magnetic field, additional types of fermionic order are triggered by the condensation of the 
SVDW and CDW order parameters. To show that, we first introduce the Zeeman coupling between the uniform field 
H and the electrons: 


Mzeeman — ^ H 


k.z 


(B16) 


We also introduce the Qi + Q 2 = (tt, tt) charge-current density-wave A^ (i.e. a purely imaginary CDW) and the 
spin density-wave A'g (i.e. a purely real SDW) defined by: 
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'^S ^ ^ ' ^aj3 ^^62 ,k:a*^ei ,k^ *^ei ,ka*^e2 ,k^^ 

k 

1-iiC = Ac'^a/3 ('^l2,ka^ei,k/3 ~ ,ka^e2,k/3) 

k 

Following the same steps as in the previous subsection, we obtain the expanded action: 

^eff [Ml, M2, A^ 5 , A;i] ^ ^ ^Tr [^0 (Vmag + ^5 + V^C + V^Zeeman) 

n 

where the Nambu-space matrices are given by: 


Vze 


and: 


as well as: 



/ 

-H, 

“1“ 

0 


0 

0 0 



-H,- 

:Hy H 

z 

0 


0 

0 0 



0 

0 

-H, 


Glx T 

0 0 



0 

0 


-H,- 

iHy 

hz 

0 0 



0 

0 

0 


0 

—Hz — Hx H“ iH^ 


V 

0 

0 


0 


0 

~Hx ^Hy Hz 





/O 0 

0 

0 

0 0 \ 






0 0 

0 

0 

0 0 





V)c = 

0 0 

0 

0 

zA" 0 






0 0 

0 

0 

0 iA'l; 






0 0 

-^A'^ 

0 

0 0 






\0 0 

0 

-zA'S 

0 0 / 


/o 

0 


0 


0 


0 

0 

0 

0 


0 


0 


0 

0 

0 

0 


0 


0 


-A),.. 

“ (As,x “ *a;^ ,^) 

0 

0 


0 


0 


{^'s,x + ^A's^y) Ag^^ 

0 

0 


A's.. 

-(A's 

X - *A 

kv) 

0 

0 

VO 

0 

-(Ak 

X + ^^'s,y) 

A's.. 


0 

0 


Vs = 


A straightforward evaluation yields, to leading order in the magnetic field: 


5 eff = 5 eff [H = 0 ] + C (H-M 


(H-Ma)' 


+ 7 [A" H • (Ml X M2) + (H • A's) (Ml • M2)] 

+ 7 [(Ml • H) (M2 • A's) + (M2 • H) (Ml • A's)] 

where we neglected all isotropic biquadratic terms of the form H^Mf. The coefficients are given by: 


^2 r^2 


7 = 4 J Gh,kGei,kGe2,k {Gei,k + Ge2,k — Gh,k) 

= 4 f fcGei,feGe2,fc 
Jk 


(B 17 ) 


(B 18 ) 


(B 19 ) 


(B 20 ) 


(B 21 ) 


(B 22 ) 


r] 


(B 23 ) 
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It is useful to perform an expansion around perfect nesting, i 5 o = ^2 = 0 . The coefficients C and rj become identical 
in this limit: 


, PF (7(m 


(B 24 ) 


The fact that C > 0 implies that the magnetic field induces an easy plane, rather than an easy axis anisotropy. As 
for the coefficient 77, it remains zero in all orders in perturbation theory if an infinite bandwidth is assumed. However, 
keeping the top of the hole pocket W (or bottom of the electron pocket) throughout the calculation gives: 


PF(wy^ 


(B 25 ) 


The fact that 7 0 implies that, in the presence of a uniform field, the SVDW order parameter <PsvDW ^1 ^ ^2 

also triggers a charge-current density-wave A^, whereas the CDW order parameter :/Jcdw oc Mi • M2 triggers a spin 
density-wave of same period, Ag. Although this was expected by symmetry, here we have microscopic expressions 
for the corresponding Ginzburg-Landau coefficients. It is interesting then to compare the coefficient 7 in Eq. (IB 22 I) . 
which determines the amplitudes of and Ag, to the coefficient A in Eq. (|B 13 ( 1 . which determines the amplitudes 
of Ag and A^. We hnd that: 


jH 


- 2.3 



(B 26 ) 


Therefore, for pnictides whose band dispersions do not deviate strongly from perfect nesting, and whose bandwidths 
are not too large either, it is conceivable that the two coupling constants jH and A will be of similar order for moderate 
values of the magnetic field H. As a result, the charge-current density-wave and the spin density-wave generated in 
the presence of the field could be as large as the spin-current density-wave and the charge density-wave generated in 
the absence of the field. 


Appendix C: Superconducting pairing interactions 

Here we show explicitly that fluctuations associated with an imaginary SDW instability or with a real CDW 
instability give rise to attractive pairing interactions. For our purposes, it is sufficient to consider only the two Fermi 
pockets connected by the momentum transfer Q = (tt, tt) associated with these two ordered states. To simplify the 
notation, here we will denote the fermionic operators associated with these bands by dkcr and /ko-. In both cases, k 
is measured relative to the center of each Fermi pocket. Consider first the action describing the coupling between the 
electrons and the complex SDW bosonic held A5 = A5Z (for simplicity, we consider it polarized along the z axis): 


S = - 


I (*w„ £^d,k) <^ko-<^k(T + (*w„ £/,k) /kcr/kcr + / O' ( Ag -|- i^s,-k-qfka^' 

Jk Jk,q '• 

+ [ Xs^ + f x-g (k,H„) 

Jk Jk 


^S,k^S,-k 


(Cl) 


where k = (k, ojn), = T f ^ (with the appropriate bosonic fin or fermionic a;„ Matsubara frequency), and we 

n ' ' 

left implicit the sum over spin indices, as well as the dependence of the fermionic operators on the fermionic Matsubara 
frequencies, xs and XiS are the susceptibilities associated with the real and imaginary SDW, and g is the coupling 
constant. Note that, because Q is a commensurate vector, the real and imaginary SDW helds are independent. 
Introducing the four-dimensional Nambu operator: 


^k ( '^kt /k-l" /-kt 'j 


(C 2 ) 
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the action can be written conveniently as: 

S = — (iWni — £k) 'I'k + f Xs + [ XiS 

Jk Jk Jk 

~^9 f ^'s-k-q'^i.PS'^p + 9 f ^S-k-q'^kPiS'^p 

k.Q J k,Q 


A" A" 
^S,k^S-k 


k,q 

where we defined the 4x4 matrices: 


(C3) 


£k = ( ^ j ; PS = ( ^ ^° I = To 0 cr^ ; PiS = ( *5* = -r^ 0 CTy 

V 0 EfMTz I \ To 0 / \ -ITz 0 / 


(C4) 


where Ti are the Pauli matrices and 0 denotes the 2x2 matrix whose elements are all zero. To obtain the Eliashberg-like 
gap equations, we need to solve Dyson’s equation: 




with Gq J, = iujnl — Ek S'lid the one-loop self-energy: 


E. = 


9^ / Xs{k- q) PsGqPs + 5^ / XiS (k - q) PisGqPis 

J a J a 


It is convenient to parametrize the self-energy by: 


~ ~ a) ” 1 ” As + ^k 


(C5) 


(C 6 ) 


(C7) 


where we introduced the imaginary normal components the real normal components and the anomalous 

components Wfj^^k (p = d, / is a band index): 


Zk = 


Zd,kTo 0 
0 Zf^kTo 


Wk = 


Wd,kTx 0 

0 Wf,kTx 


id,kTz 0 


(C 8 ) 


The superconducting gap in band p is therefore proportional to Wf^^k- Using Eqs. (IC5II and (IC7E it is straightforward 
to invert the matrix and obtain G. Substituting it in (IC 6 I) and comparing back with Eq. (IC7E we arrive at a set of 
six self-consistent equations. Four of them have the same form for either real or imaginary SDW, namely, the two 
equations that renormalize the dispersion ^ = ^a,k + Sk and the two that renormalize the quasi-particle weights 
Z^.^-- However, the two self-consistent gap equations acquire different forms: 


Wd,k = - [ [ff As {k - ?)] - [ [-9^X^S {k - ?)] 

jq ^f,q Jq 

[ff As {k - g)] ^ [-pAiS (fc - g)] 


WfM = - 


Df,q 

Wd,q 

Dd,q 


(C9) 


where we defined I’'™™ form of these equations, it becomes clear that while the 

fluctuations near the real SDW instability give rise to a repulsive inter-band pairing interaction, Vdf oc q^xSj the 
fluctuations near the imaginary SDW instability promote an attractive inter-band pairing interaction, Vdf oc —g^XiS- 
This difference relies ultimately on the different structures of the matrix elements (IC4I) in Nambu space. 

A similar analysis can be performed in the charge channel: 


S = - 


/ {i<-kln Ed.k) dAdko- + {i^n E/^k) /A/ko- 9 (_^_gdA/q(T + ^C,-k-qfka-^<l<^) 

Jk Jk,q ^ ^ 

+ [ Xc^ (k, Xln) Ap fcAp -I- f Xic (h; H„) ■ 

Jk Jk 


I ^C,k^C,-k 


(CIO) 
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In Nambu space, we obtain: 


S — — f (jWni — Ek) 'I'k + f Xc^ ^c,k^'c,-k + [ XiC 0^: ^n) ^C^k^C-k 

J k J k J k 

+9 ( ^'c,-k-q'^i.PC'^p 9 f ^C,-k-q'^kP^C'^p 

J k,q J k,q 

where we defined the 4x4 matrices: 

f 0 Tz \ _ ^ - f 0 iro \ „ 

Pc = \ ^ ] = ; p,c = [ ^ ]= -To ® Cy 

\Tz 0 J y -iTo 0 J 

Solving the one-loop Dyson equation, we obtain the two self-consistent gap equations: 


Wd,k 

Wf,k 



Therefore, in the charge channel, real CDW fluctuations promote inter-band pairing attraction, Vdf oc 
whereas imaginary CDW fluctuations promote repulsion, Vdf oc g^XiC- 


(Cll) 

(C12) 

(C13) 

-g^xc, 







